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ABSTRACT 

A-optimal  and  MV  -optimal  repeated  measurements  designs  are  given  both  for  direct 
and  residual  treatment  effects,  for  comparing  several  test  treatments  with  a  control.  The 
models  considered  are  basically  of  two  types:  without  preperiods  and  the  circular  model.  It 
is  shown  that  some  known  balanced  and  strongly  balanced  uniform  repeated  measurements 
designs  can  be  modified  to  obtain  optimal  designs  for  this  problem.  Some  other  methods 
of  finding  optimal  designs  are  also  given. 

1.  INTRODUCTION 

We  consider  the  problem  of  finding  optimal  experimental  designs  for  comparing  t  test 
treatments  with  a  control,  using  n  experimental  units  which  are  observed  in  each  of  p 
periods;  the  model  is  the  homoscedastic,  additive,  repeated  measurements  linear  model. 
Recently,  this  problem  has  been  considered  by  Pigeon  (1984)  and  Pigeon  and  Raghavarao 
(1987).  They  have  defined  classes  of  ‘good’  designs  and  investigated  their  efficiencies. 
Though  they  did  not  consider  the  problem  of  finding  optimal  designs,  we  shall  show  that 
some  of  their  designs  are  indeed  optimal. 

The  problem  of  determining  optimal  designs  in  repeated  measurements  models  where 
all  the  elementary  treatment  contrasts  are  equally  important,  i.e.,  there  is  no  special  treat¬ 
ment  like  a  control,  has  been  investigated  by  several  authors.  Recent  research  includes 


Hedayat  and  Afsarinejad  (1975,  1978),  Cheng  and  Wu  (1980),  Magda  (1980),  Kunert 
(1983,  1984a,  1984b,  1985),  Mukhopadhyay  and  Saha  (1983),  Sen  and  Mukherjee  (1987), 
Afsarinejad  (1983,  1985)  and  Constantine  and  Hedayat  (1982).  Hedayat  and  Afsarinejad 
(1975)  has  an  excellent  bibliography  of  eariler  literature. 

The  existing  literature  on  optimal  designs  for  comparing  treatments  with  a  control 
deal  only  with  additive,  homoscedastic  models  eliminating  0, 1,  or  2  ways  of  heterogeneity, 
with  no  residual  effects.  Even  though  this  literature  does  not  have  any  optimal  design 
for  comparing  treatments  with  a  control  in  repeated  measurements  models,  some  methods 
which  have  been  used  there  will  prove  to  be  valuable  in  the  present  context.  One  of 
the  relevant  methods  is  due  to  Notz  (1985).  He  showed  that  if  one  starts  with  a  latin 
square  of  order  t 2  + 1,  using  the  symbols  1,2, . . . ,  t2  + 1,  and  redefines  each  of  the  symbols 
t2  +  1, . . . ,  t2  +  t  as  0,  which  denotes  the  control,  and  calls  the  symbols  1, . . .  t2  the  test 
treatments,  then  the  resulting  design  is  optimal  for  comparing  t 2  test  treatments  with  a 
control  in  additive  models  eliminating  heterogeneity  along  t2  + 1  rows  and  t2  + 1  columns. 
It  is  interesting  to  note  that  the  initial  latin  square  was  itself  universally  optimal  (Kiefer, 
1975)  for  comparing  t2  +  t  treatments  among  themselves;  all  treatments  being  equally 
important  in  the  criterion. 

In  section  3  of  this  paper  we  successfully  use  this  type  of  a  technique  in  the  more  com¬ 
plicated  setup  of  repeated  measurements  models.  More  precisely,  we  show  that,  in  these 
models,  an  optimal  design  for  comparing  several  treatments  with  a  control  can  sometimes 
be  obtained  by  starting  with  a  suitable  optimal  design  for  the  traditional  setup  of  compar¬ 
ing  all  treatments  among  themselves  (with  no  special  treatment  like  a  control)  and  then 
redefining  some  treatment  symbols  to  be  the  control.  Many  optimal  designs  to  be  given 
in  section  3  will  be  of  this  type.  It  is  interesting  to  note  that  Pigeon  (1984)  had  suggested 
this  as  one  of  the  methods  of  getting  ‘efficient’  designs. 

Our  method  of  proof  is  different  from  that  of  Notz  (1985).  We  shall  essentially  cou¬ 
ple  some  orthogonality  type  conditions  with  optimality  in  simpler  models  to  prove  our 
results.  This  technique  had  been  used  by  Kiefer  (1975)  for  proving  the  optimality  of  reg¬ 
ular  generalized  Youden  designs,  by  Magda  (1980),  Kunert  (1983,  1984a,  1984b)  and  Sen 
and  Mukherjee  (1986)  for  repeated  measurements  models  and  by  Jacroux  (1987b)  and 


Majumdar  (1986)  for  comparing  treatments  with  controls  in  additive  classification  models 
eliminating  heterogeneity  in  two  directions.  In  the  repeated  measurements  setup,  the  most 
extensive  use  of  this  technique  has  been  made  by  Kunert  in  the  papers  mentioned  above. 

We  devote  section  2  to  preliminary  notations,  definitions  and  basic  lemmas.  The 
optimal  designs  for  repeated  measurements  models  are  given  in  sections  3  and  4.  Section 
3  contains  optimal  designs  which  can  be  obtained  by  the  technqiue  due  to  Notz,  while 
section  4  briefly  outlines  some  other  techniques.  Our  results  are  illustrated  by  examples. 

2.  PRELIMINARIES 

We  are  interested  in  an  experiment  in  which  n  individuals  or  units  are  to  be  observed 
in  each  of  p  periods  upon  application  of  one  of  t  treatments  in  each  period.  An  allocation 
of  the  treatments  is  called  a  repeated  measurements  design,  to  be  denoted  by  RMD(t,n,p). 
The  set  of  all  RMD(t,n,p)  will  be  denoted  by  0(t,n,p).  If  d  is  an  RMD(t,n,p),  then  d(i,j) 
will  denote  the  treatment  assigned  to  period  i  of  unit  j.  The  corresponding  observation 
will  be  denoted  by  Yu, . 

Suppose  t  test  treatments  are  to  be  compared  with  a  control,  which  is  a  standard 
treatment,  and  a  design  d  is  selected  from  VI  (t  +  l,n,p).  Then  the  model  is 


Ydi)  =  P  +  <*%+  0]  +  rd(iJ)  +  Pd(i~lJ)  +  *»;•  (2.1) 

Here  *  =  1, . . . ,  p,  j  =  1, . . . ,  n;  d{i,j)  =  0, 1, . . . ,  t  with  0  denoting  the  control  and  1, . . . ,  t 
denoting  the  test  treatments.  The  usual  assumptions  regarding  e  holds  in  this  model,  i.e., 
E(e)  =  0,  V(c)  =  o2,  e’s  are  uncorrelated.  The  symbols  p  stands  for  a  general  effect,  a 
for  the  period  effect,  0  for  the  unit  effect,  r  for  the  direct  effect  of  the  treatment  and  p  for 
the  residual  effect  of  the  treatment  applied  in  the  previous  period.  When  Pd(oj)  =  0  for 
all  j,  then  this  the  is  model  without  preperiods  (see  Hedayat  and  Afsarinejad  (1975,  1978), 
Cheng  and  Wu  (1980)).  On  the  other  hand  it  is  sometimes  convenient  to  have  a  period 
preceeding  period  1.  The  observations  in  this  period  are  not  used  in  the  analysis;  the 
treatments  are  applied  principally  to  generate  residual  effects  in  period  1.  Among  models 
with  preperiods,  we  shall  consider  only  the  circular  model  proposed  by  Magda  (1980).  Here 
=  PH,j)  for  all  j. 


If  Yd  denotes  the  vector  of  all  np  observations,  then  as  in  Kunert  (1983),  we  may  write 

Yd  =  InpM  +  Pa  +  Ufi  +  Tdr  +  Fdp  +  c  (2.2) 

where  la  denotes  an  a  x  1  vector  of  unities,  a  =  (at,. . .  ,  ap)',  /?  =  (fiu. . .  ,/?„)',  r  = 
(r0, . . . , rty,  p  ■■  (po. •  •  ■ , PtY .  €  =  (ci, . . . ,  Crip)'".  U,  Td  and  Fd  having  obvious  interpre¬ 
tations.  Let  us  define, 

ftrftu  =  number  of  appearances  of  treatments  i  on  unit  u 

ndiu  =  number  of  appearances  of  treatment  t  on  unit  u  in  the  first  p  —  1  periods 

f-dik  =  number  of  appearances  of  treatment  i  in  period  k 

mdij  —  number  of  appearances  of  treatment  i  preceeded  by  treatment  j  in  the  same  unit 

rdi  =  number  of  appearances  of  treatment  t 

fdi  —  number  of  appearances  of  treatment  i  in  the  first  p  -  1  periods. 

Before  proceeding  to  find  optimal  designs  for  comparing  treatments  with  a  control, 
we  have  to  define  some  designs  which  are  optimal  in  the  setup  where  all  treatments  are 
equally  important,  without  any  special  treatment  like  a  control.  To  do  this,  we  briefly  (for 
Definitions  2.1,  2.2  and  2.3  only)  depart  from  our  established  notation  and  number  the 
treatments  in  Cl(t,n,p)  by  1, . . . ,  t;  0  is  absent  here  since  there  is  no  control. 

DEFINITION  2.1.  (Heda.ya,t  and  Afsarinejad  (1975)).  An  RMD(tji,p)  is  a  balanced  uni¬ 
form  repeated  measurements  design,  denoted  by  BURMD(t,n,p ),  if 


Tdl  ^dt 

(2.3) 

ridu  =  ...  =  ndin  for  each  i, 

(2.4) 

fdii  ——  ...  =  tdip  for  each  t, 

(2.5) 

and 


mdij  =  md\7  for  all  i,j,i  ^  ;. 


DEFINITION  2.2.  (Cheng  and  Wu  (1980)).  An  RMD(t,n,p)  is  a  strongly  balanced  uniform 
repeated  measurements  design,  denoted  by  SBURMD(t,n,p),  if  (2.3),  (2.4),  (2.5)  hold  and 

mdij  =  mdn  for  ail  i,j. 

DEFINITION  2.3.  (Magda  (1980)).  An  RMD(tji,p)  is  a  circular  strongly  balanced  uni¬ 
form  repeated  mesurements  design,  denoted  by  CSBURMD(tji,p),  if  (2.3),  (2.4),  (2.5) 
hold  and 

tdij  =  «<m,  for  all  i,j, 

where 

edij  =  TTldij  +  fdij 

with 

fdij  =  number  of  units  u  with  d(l,u )  =  i  and  d(p,  u)  =  j. 

There  are  several  results  establishing  the  universal  optimality  for  direct  as  well  as 
residual  treatments  effects  of  these  RMD's.  Hedayat  and  Afsarinejad  (1978)  proved  that 
a  BURMD  is  universally  optimal  in  the  class  of  all  designs  with  the  qualification  that 
no  treatment  appears  in  consecutive  periods.  Cheng  and  Wu  (1980)  proved  that  a  design 
obtained  by  repeating  the  observations  in  the  last  period  of  a  certain  BURMD  is  universally 
optimal.  Cheng  and  Wu  (1980)  also  proved  the  universal  optimality  of  an  SBURMD.  All 
these  results  are  for  a  model  without  preperiods.  On  the  other  hand,  Magda  (1980)  proved 
the  universal  optimality  of  a  CSBURMD  in  a  circular  model.  There  are  many  more  results 
on  optimality  of  RMD’s  in  these  and  the  other  papers  mentioned  in  the  introduction. 

In  case  all  the  treatments  are  equally  important,  the  parametric  functions  of  interest 
are  usually  an  orthonormal  basis  of  the  treatment  contrasts.  On  the  other  hand,  in  the 
case  where  several  test  treatments  (l,...,t)  are  being  compared  with  a  control  ( o ),  the 
parametric  functions  of  direct  treatment  effects  which  are  of  interest  are 


(r0  —  ri), . . . ,  (to  -  rt).  These  are  estimated  by  their  BLUE’s  (tdo  —  ^ai), . . . ,  {?do  -  fat)* 
under  the  design  d.  The  problem  here  is  to  choose  a  design  for  which  these  estimates 
are  “most  accurate”.  Sometimes  the  object  of  an  experiment  may  be  to  estimate  the 
corresponding  residual  treatment  contrasts,  (p0  —p j), . . . ,  (po~Pt),  which  are  estimated  by 
(pdo  —  Pdi),  •  •  • ,  (Pdo  ~  Pdt)-  We  suspect  that  in  most  experiments,  these  residual  treatment 
contrasts  have  less  importance  than  the  direct  treatment  contrasts. 

Not  all  individual  criteria  belonging  to  the  class  which  forms  universal  optimality 
may  be  statistically  meaningful  in  the  context  of  comparing  test  treatments  with  a  con¬ 
trol.  There  are,  however,  two  criteria,  A-optimality  and  MV  -optimality,  which  do  possess 
natural  statistical  interpretations,  and  have  received  considerable  attention  in  the  case  of 
optimal  block  designs  and  optimal  row-column  designs.  These  are  given  in  Definitions  2.4 
and  2.5. 

DEFINITION  2.4.  (i)  do  E  n(t  +  l,n,p)  is  A-optimal  for  the  direct  treatment  effects  if 

t  t 

£  Var(f*,o  -  f*,i)  <  £  Vartfjo  -  h.)  (2.6) 

»=1  1=1 

for  all  d  E  fi(t  +  I,n,p). 

(ii)  do  €  fl (t  +  l,n,p)  is  A-optimal  for  the  residual  effects  of  the  treatments  if  (2.6)  hold 
for  all  d  £  fl(i  +  l,n,p)  with  f’s  replaced  by  p*s. 

DEFINITION  2.5.  (i)  do  e  n(t  +  l,n,p)  is  MV-optimal  for  the  direct  treatment  effects  if 


Max{  Var(tdoO  ~  fd0i)  ■  *  =  1, . . . ,  t)  <  Max{  Var(fdo  —  ?*, :  *  =  1, . . . ,  t) 


(2.7) 


for  all  d  G  n(t  +  l,n,p). 

(ii)  do  E  Q (t  +  l,n,p)  is  MV -optimal  for  the  residual  effects  of  the  treatments  if  (2. 7)  holds 
for  all  d  E  0(t  +  l,n,p)  with  f’s  replaced  by  p’s. 

The  purpose  of  this  paper  is  to  find  A-  and  MV-optimal  designs.  We  shall  start  with 
two  basic  lemmas.  Consider  a  class  of  designs  P  and  two  Gauss-Markov  models  for  a 
dEP ; 

6 


E(Yd)  =  XldOt  +  X2d02  +  X3d03  V{Yd)  =  a2 1 


(2.8) 


E(Yd)=Xld01  +  X2d0  2  V(Yd)=o2I  (2.9) 

The  model  (2.8)  has  been  called  ‘finer’  than  (2.9)  by  Kunert  (1983). 

Our  parameter  vector  of  interest  is  0\;  more  precisely,  we  would  like  to  estimate 
the  vector  QO i  where  Q  is  a  matrix.  The  vectors  02,03  and  the  scalar  a2  are  nuisance 
parameters.  Let 


Ad  =  X'ld(I  -  X2d(X'uX2d)-X'2i)Xu 


and 


cd  =  x[d{i  -  zd{z,dzdyz,d)xld 


where 


Zd  —  {X2d  :  ^3(i)- 


Suppose  is  a  real  valued  function  defined  on  nonnegative  definite  matrices.  Then  do  €  P 
is  V'optimal  for  estimating  Q6\  under  (2.8)  if 


<  rJ>{QCjQ')  for  all  deP, 

while  it  is  V’-optimal  for  estimating  QO i  under  (2.9)  if 


'l>{.QAdaQ')  ^  'l>{QAd  Q')  for  a11  dep. 


To  avoid  trivialities,  we  shall  assume  that  Q0\  is  estimable  for  each  d  G  P  under  (2.8). 
Hence  Q0\  is  estimable  under  (2.9)  also.  We  also  assume  that 


Cd  1  =  Adl  —  0  for  all  d  €  P, 


where  1  is  a  vector  of  unities.  Consequently  <J1  =  0. 

The  question  is:  when  is  a  design  which  is  0-optimal  under  (2.9),  0-optimal  under 
(2.8)  also?  An  answer  is  given  in  Lemma  2.1. 

LEMMA  2.1.  Let  0  have  the  property  that  0(A  —  B)  <  0(A)  whenever  A,  B  and  A  -  B 
are  nonnegative  deBnite.  Suppose  do  6  P  satisBes 

x'^x^  =  X'^X^X'^X^rx'^X^  (2.10) 

and 

HQA70Q')  <  0(GAjg')  for  alldep  (2.11) 

then  d0  is  0- optimal  for  estimating  Q0\  under  (2.8). 

PROOF:  By  Lemma  3.1  of  Magda  (1980)  (also  reproduced  in  Proposition  2.3  of  Kunert 
(1983)),  (2.10)  implies 


Cdo  —  Ado  ■ 

Hence 

'KQCToQ')  =  WQAToQ') 

<  0(QAJQ')  by  (2.11) 

<  0(QCd-Q')  for  all  d  6  V 

since  there  exists  ^-inverses  such  that  Cj  -  Aj  is  nonnegative  definite,  because  Ad  -  Cj. 
is  nonnegative  definite  (cf.  Wu  (1980)).  Hence  the  lemma. 

For  the  special  case  when  X^d  =  0  for  all  d  €  D  (that  is,  8 1  is  absent  in  (2.8)  and 
(2.9)),  and  X\dX\d  is  nonsingular  for  all  d  €  P,  a  somewhat  stronger  result  is  available 
in  Theorem  3.1  of  Majumdar  (1986).  We  reproduce  this  as  Lemma  2.2,  for  the  sake  of 
completeness. 


LEMMA  2.2.  Let  0  have  the  property  that  0(A  —  B)  <  0(A)  whenever  A,  B  and  A  —  B 
axe  nonnegative  deSnite  matrices.  Suppose  do  G  V  satisfies 


*k,xuo(xu*M.)"‘e'  =  ° 


(2.12) 


and 


HQ(X[doXldo)-lQ,)<tj,(Q(X,ldXld)-1Q,)t  for  all  deD. 
Then  do  is  0- optimal  for  estimating  QBi  under  (2.8). 


(2.13) 


Examples  of  functions  0  which  satisfy  the  conditions  of  Lemmas  2.1  and  2.2  are: 

0  =  trace,  giving  the  A-optimality  criterion  and  0  =  maximum  diagonal  element,  giving 
the  M K -optimality  criterion. 


3.  OPTIMAL  DESIGNS 


In  this  section  we  first  apply  Lemma  2.2  to  the  repeated  measurements  model  (2.1) 
without  preperiods.  Here 


D  =  n(t  +  l,n,p) 


Q  = 


( 1-1  0 
1  0  -1 


0  \ 
0 


Vi  o  o  ...  -lj 

atxt+1  matrix;  and  equating  models  (2.2)  and  (2.8)  we  substitute, 


Xid  =  T4,  Xzd  =  (lnp  :  P  :  U  :  Fj]. 


Recall  that  in  our  notation  Xjd  =  0  for  all  d  €  V.  To  find  A-optimal  designs,  we  choose 
0  =  trace.  Hence,  condition  (2.13)  is  satisfied  if 


i  *  ^do  1 1  r*,o  rdo  l  ^/t. 


(3.1) 
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Condition  (2.14)  tells  us  that  if  we  denote 


H  =  . 

then 

1  npTdoH  =  0,  P'T^H  =  0,  U'T^H  =  0,  F^H  =  0. 
These  conditions  lead  us  to  the  following  relations: 


^  do  lk 

.  •  •  f-dotky  f-ioOk  t-dq  lky/t  for  k  —  1,  •  •  •  ,  P 

(3.2) 

Tldolu  = 

. . .  ftdoOu  ~~  ftdo xuVt  for  u  =  l,...,n 

(3.3) 

mdolO  ~ 

. . .  =  m^toyin^oo  =  m^io-^i 

(3.4) 

mdo  It 

. . .  —  tn\d0ti\  rndoOi  =  y/t  for  J  —  1, . .  • ,  t. 

(3.5) 

This  shows  that  any  RMD(t+l,n,p)  satisfying  relations  (3.1)  -  (3.5)  is  ^-optimal  for  direct 
treatment  effects  in  H(t  +  l,n,p)  for  a  model  without  preperiods.  These  designs  are  MV  - 
optimal  also,  since  (3.1)  satisfies  (2.13)  with  tp  =  maximum  diagonal  element.  One  way  to 
obtain  these  optimal  designs  is  given  in  Theorem  3.1  which  can  be  established  by  verifying 
conditions  (3.1)  -  (3.5). 

THEOREM  3.1.  Let  d*  be  a  strongly  balanced  uniform  repeated  measurements  design 
SBURMD  (to2  +  u;,n,p),  where  w  is  a  positive  integer.  Let  do  be  obtained  from  d*  by 
changing  each  of  the  treatment  symbols  w2+l,...,w2  +  w  to  the  symbol  0  which  denotes 
the  control ,  and  keeping  everything  else  unchanged.  Then  do  is  A-  and  MV -optimal  for 
direct  treatment  effects  for  comparing  w 2  test  treatments  with  a  control  in  D(w2  +  l,n,p) 
for  a  model  without  preperiods. 

Example  3.1.  Methods  for  constructing  SBURMD’s  have  been  given  in  Cheng  and  Wu 


i 

k 

U 

i 


(1980).  Using  these,  we  give  an  A-  and  MV-optimal  design  for  direct  treatment  effects  in 
0(5,36, 12)  in  Table  1.  Here  tv  =  2,  n  =  36,  p  =  12;  0  denotes  the  control  and  1,2, 3, 4  the 


test  treatments. 


MS 


Table  1:  Showing  the  optimal  design  of  Example  3.1 


Periods 


E 


ft 


Units 


1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 
4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


1 

2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 


2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 
4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
1 
1 
1 
1 
1 
1 


2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 


3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 
4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
2 


3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 


4 

4 

4 

4 

4 

4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 


4 

0 

0 

1 

2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 
3 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 
4 


0 

0 

1 

2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 


0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

4 
4 
4 
4 
4 
4 
0 
0 
0 
0 
0 
0 


0 

1 

2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 
0 
1 
2 

3 

4 
0 


11 


Let  us  now  turn  our  attention  to  the  estimation  of  residual  treatment  effects  in  models 
without  preperiods.  Let  us  substitute  in  Lemma  2.2, 

Xid  =  fdt  Xsd  =  [lnp  •  P  •  U  :  Td\. 

Condition  (2.13)  is  satisfied  for  both  A-  and  AfV -optimality  if 

'doi  •••  —  y/t  (3.6) 

Then  condition  (2.12)  gives  us  the  following  relations 


lk  •  •  •  —  t-dotki  t-doOk  —  ^  Ik  ^  for  k  —  1, . .  .  ,  p  1  (3-7) 

=  •••  =  ^(mAdoOii  "  ri^oiM  yft  for  u  =  1, . . .  ,n  (3.8) 

m  da  01  =  ...  =  rrido01tmdo  00  —  mde01  >ft  (3.9) 

“  ...  rn^oto  y/t,  i  —  1, . . . , t  (3.10) 


Hence  the  conditions  (3.6)  -  (3.10)  are  sufficient  for  a  design  do  to  be  A-  and  MV- 
optimal  for  residual  effects  in  fl(t  +  l,n,p).  The  design  do  of  Theorem  3.1  does  not  satisfy 
(3.8).  A  family  of  optimal  designs  which  do  satisfy  these  conditions  are  given  in  Theorem 
3.2. 

THEOREM  3.2.  Let  d*  be  obtained  by  repeating  the  observations  in  the  last  period  of 
a  balanced  uniform  repeated  measurements  design  BURMDfw 3  +  w, X(w2  +  w),w2  +  u>), 
where  A  and  w  are  two  arbitrary  positive  integers.  Let  do  be  obtained  from  d*  by  changing 
each  of  the  treatment  symbols  w2  +  1, . . . ,  w3  +w  to  the  symbol  o  which  denotes  the  control 
and  keeping  everything  else  unchanged.  Then  do  is  A -  and  MV -optimal  for  residual 
treatment  effects  for  comparing  w 2  test  treatments  with  a  control  in  Cl(w2  +  l,A(w3  + 
u>),u/3  +  w  +  1)  for  a  model  with  preperiods. 


To  verify  conditions  (3.6)  -  (3.10)  it  helps  to  realize  that  in  the  original  BURMD(u>2  + 
w,  A(u>2  +  w),w7  4-  to)  each  treatment  occurs  X  times  in  each  period  and  each  treatment 
is  preceeded  by  every  other  treatment  X  times  (see  Hedayat  and  Afsarinejad  (1975,  p. 
231).  Cheng  and  Wu  (1980, Corollary  3.3.1)  established  optimality  properties  of  a  BURMD 
augmented  by  it  last  period,  for  comparing  all  treatments  among  themselves,  no  treatment 
being  a  control. 

Example  2.3.  Starting  with  the  BURMD(6,6,6)  from  Hedayat  and  Afsarinejad  (1975, 
p.  232)  we  give  an  A-  and  MV- optimal  design  for  residual  treatment  effects  in  fi(5,6,7). 
Here  the  control  is  denoted  by  0  and  the  4  test  treatments  by  1,  2,  3  and  4;  w  =  2  and 
X  =  1. 

Units 

1  2  3  4  0  0 

0  1  2  3  4  0 

2  3  4  0  0  1 

Periods  0  0  1  2  3  4 

3  4  0  0  1  2 

4  0  0  1  2  3 

4  0  0  1  2  3 

Finally,  let  us  consider  a  repeated  measurements  model  which  is  circular.  For  direct 
treatment  effects  we  again  start  with  the  substitution  Xu  —  Tj,  and  Xu  =  [lnp  :  P  : 
U  :  F,*].  It  is  clear  that  a  sufficient  condition  for  do  to  be  A-  and  MV  -optimal  is  that 
the  conditions  (3.1)  -  (3.5)  are  satisfied,  with  the  exception  that  md%j  is  replaced  by  e*/. 
The  symbol  Cdij  has  been  defined  in  Definition  2.3  -  it  is  the  number  of  times  treatment 
*  follows  treatment  j,  counting  in  a  circular  fashion.  On  the  other  hand,  we  substitute 
X\ d  =  Fd  and  Xu  =  (lnp  :  P  :  U  :  Td\  to  get  optimal  designs  for  residual  treatment 
effects.  The  sufficient  conditions  for  the  optimality  of  do  are  (3.1),  (3.2),  (3.3)  and  (3.9), 
(3.10)  with  mdij'a  replaced  by  s.  Families  of  optimal  designs  are  given  in  the  following 
theorem,  which  is  easily  established  by  verifying  the  sufficient  conditions  mentioned  above. 

THEOREM  3.3.  Let  d*  be  a  circular  strongly  balanced  uniform  repeated  measurements 
design  CSBURMDftv 2  +  w,n,p),  where  w  is  a  positive  integer.  Let  do  be  obtained  from 


d*  by  changing  each  of  the  treatment  symbols  tv2  +  1, . . . ,  w2  +  w  to  the  symbol  o  which 
denotes  the  control,  and  keeping  everything  else  unchanged.  Then  do  is  A-  and  MV - 
optimal  for  direct  as  well  as  residual  treatment  effects  for  comparing  w2  test  treatments 
with  a  control  in  fl(w3  +  l,n,p)  for  a  circular  model. 

Example  3.3.  Using  the  methods  for  constructing  CSBURMD’s  given  by  Sen  and  Mukher- 
jee  (1985),  we  give  an  A-  and  M V -optimal  design  in  0(5,6, 12),  under  the  circular  model 
for  direct  as  well  as  residual  treatment  effects.  Here  w  =  2;  0  denotes  the  control,  while 
1,2, 3, 4  denote  the  4  test  treatments 


Periods 


1 

0 

2 

0 

3 

4 

4 

3 

0 

2 

0 

1 

2 

1 

3 

0 

4 

0 

0 

4 

0 

3 

1 

2 

Units  3 

2 

4 

1 

0 

0 

0 

0 

1 

4 

2 

3 

4 

3 

0 

2 

0 

1 

1 

0 

2 

0 

3 

4 

0 

4 

0 

3 

1 

2 

2 

1 

3 

0 

4 

0 

0 

0 

1 

4 

2 

3 

3 

2 

4 

1 

0 

0 

4.  SOME 

MORE 

OPTIMAL 

DESIGNS 

In  section  3  we  applied  Lemma  2.2  to  obtain  optimal  designs.  In  this  section  we  shall 
apply  Lemma  2.1.  The  optimal  designs  in  this  section  cannot  be  obtained  by  Notz-type 
methods,  in  general. 

The  first  step  consists  in  equating  models  (2.2)  and  (2.8).  This  means  that  we  have 
to  partition  the  set  of  matrices 


{ifipi  P,U,Ti,Fi) 


as 


{A’lrf,  Xidi  Xsd}' 


There  are  many  possible  choices,  depending  on  how  we  partition.  We  shall  consider  only 
one  partition.  This  application  will  demonstrate  the  use  of  Lemma  2.1. 


3 

\ 

*> 

! 


Thus  we  are  looking  for  A-  and  MV -optimal  designs  for  direct  treatment  effects.  We  shall 
consider  a  circular  model  for  the  repeated  measurements. 

For  a  do  6  0(t  +  l,n,p),  condition  (2.10)  gives  us 

tidoii  =  . . .  =  ndoin  for  each  s  =  0, 1, . . .  ,t  (4.1) 

and 

edo ij  =  n  [^doil^dojp  "i”  ^doi^dojl  +  •  •  •  +  ^-doip^dojp—  l] 

for  each  t  =  0, 1, . . .  ,t  and  j  =  0, 1, . . .  ,t,  (4.2) 

where  the  notation  t  was  defined  in  Definition  2.3. 

For  tp  =  trace  or  tp  —  maximum  diagonal  element,  condition  (2.11)  says  that  do  is  A- 
or  MV -optimal  for  comparing  test  treament  with  a  control  in  the  model 

E{yik)  =  Ti  +  pk,  V (yik)  =  a2,  y,fc’s  uncorrelated. 

This  is  an  additive,  homoscedastic  model  with  treatment  and  block  effects,  the  blocks  being 
the  periods.  Many  results  giving  optimal  block  designs  for  comparing  test  treatments  with 
a  control  sure  available  in  the  literature.  To  get  started  in  this  area  the  reader  may  look  at 
Bechhofer  and  Tamhane  (1981),  Majumdar  and  Notz  (1983),  Giovagnoli  and  Wynn  (1985), 
Constantine(l983),  Hedayat  and  Majumdar  (1985),  and  Jacroux  (1987a).  This  is  a  list  of 
only  a  few,  not  all,  important  papers.  It  is  interesting  to  note  that  conditions  (2.11)  and 
(4.1)  show  that  do  is  an  optimal  design  for  comparing  test  treatments  with  a  control  in  an 
additive  model  eliminating  2  ways  of  heterogeneity  along  periods  and  units  (see  Jacroux 
(1987b,  Theorem  3.4)). 

Example  4.1.  Let  us  give  an  example  to  illustrate  this  method.  We  consider  the  class 
0(4,2,432).  There  are  3  test  treatments  (labelled  1,  2  and  3)  and  a  control  (labelled  o). 
To  display  the  design  as  an  array  let  us  define  a  2  x  24  matrix  using  the  symbols  a,  6,  c,  d: 


B(a  ^c<j\-(a^^c^^^aaeddabbdcebaadcc 
'  \b  a  a  d  e  e  a  b  b  d  c  e  b  a  a  e  d  d  a  b  b  c  d  d 

Now  consider  the  2  x  432  array  given  as 


A  =  (D,  Et  F) 


where 


D  =  (B(l,  0, 2, 0),  5(1,0, 2, 0),fl(l,0, 2, 0),£(1, 0,3,0),  B(l,  0,3,0),  5(1, 0,3,0)) 

E  =  (fl(2, 0,3,0),  B(2,0>310),B(2t0, 3,0),  B(l,  0,1,2),  B(l,  0,1,3),  J9(l,  0,2, 3)) 

F  =  (£(2,0, 1,2),  £(2,0, 1, 3),  5(2, 0,2, 3),  27(3,0, 1, 2),  5(3,0, 1,3),  i3(3,0, 2,3)). 

The  array  A,  with  rows  denoting  units  and  columns  denoting  periods  gives  an  A-  and 
MV -optimal  design  for  direct  treatment  effects  for  comparing  3  test  treatments  with  a 
control  in  0(4,2,432)  for  a  circular  model.  To  prove  the  optimality  of  A,  first  observe 
that  condition  (2.11)  for  this  setup  follows  from  Theorem  3.1  of  Hedayat  and  Majumdar 
(1984).  Finally,  the  conditions  (4.1)  and  (4.2)  can  be  established  directly.  It  can  be  seen 
that  for  this  design  do,  the  matrix 


(  (  e«*o  )  )  “ 


144  60  60  60 
60  68  26  26 
60  26  68  26 
60  26  26  68 


where  the  rows  and  columns  are  ordered  0,1 ,2,3. 

An  optimal  design  given  by  conditions  (2.11),  (4.1)  and  (4.2)  enjoy  the  model  robust¬ 
ness  property  that  they  remain  optimal  even  when  the  unit  effects  and/or  the  residual 
effects  are  zero.  The  designs  in  section  3  also  have  similar  properties.  For  example,  the 
optimal  design  do  of  Theorem  3.1  remains  optimal  when  some  or  all  of  the  following  hold: 

O: i  —  ...  Op,  ...  —  0t%  and  po  —  ...  —  . 


In  conclusion  we  observe  that  this  is  the  first  paper  which  explicitly  gives  optimal 
designs  for  comparing  test  treatments  with  a  control.  More  research  needs  to  be  done  to 
investigate  classes  fl(t,n,p)  not  covered  by  this  paper.  In  particular,  attention  needs  to 
be  directed  at  classes  where  the  number  of  periods  is  small,  since  this  is  very  useful  in 
practice.  Not  many  results  are  known  for  this  situation,  even  in  the  case  where  there  are 
no  special  treatments  like  the  control,  and  all  treatments  are  equally  important. 
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